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1. INTRODUCTION 

In the underwater engineering field, remotely operated vehicle (ROV) plays an important role in 
underwater observation, investigation, and inspection [1]—[3]. Especially in the oil and gas industry, ROV is 
used to do underwater pipe inspections as well as repairing jobs. ROV normally suffered from problems that 
include pose recovery or station keeping, underactuated conditions, coupling issues, and communication 
techniques [4]. This research paper was focusing on ROV depth control or station keeping. 

Station keeping at a certain depth is very important for underwater exploration and inspection 
[5]-[7]. However, controlling ROV is difficult because of the unexpected and unpredictable underwater 
environment [4], [8]. This is due to the nonlinear hydrodynamics effect, coupled characters of plant 
equations, lack of precise models of underwater vehicle hydrodynamics and uncertainty parameters [9], [10], 
as well as the presence of environmental disturbances [1], [11]-[14]. 

Controller design, based on simple models of underwater vehicle mass and drag, generally yields 
unacceptable performances [15]. Linear (conventional) controller is unable to adequately control the 
unmanned underwater vehicle (UUV) satisfactorily [16]. Even for a one-axis motion such as vertical motion 
or heave motion, consistent performance for a desirable range is required. Overshoot in the system is 
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unacceptable as it can harm the ROV or its inspection location [14], [17]-[20]. It is best to have as minimum 
overshoot as possible in the ROV system. 

There are many controllers designed by the researcher to cater to this problem. There are 
proportional, integral, and derivative (PID) based controllers and artificial based controllers. PID is a simple 
control technique that has been universally used because of the simplicity of implementation in a real-time 
system. Even for work class ROV, PID is used as its controller. However, the limitation is that it cannot 
dynamically compensate for unmodelled vehicle’s hydrodynamics forces or unknown disturbances. 
Furthermore, the contradictory parameter configuration such as between the rise time and overshoot may also 
exist. PID controllers have been implemented in previous work for tracking purposes in UUV [21]-[25] as 
well as in ROV [17], [26]. Normally, the PID controller was used as a basic controller to be compared with 
another complex controller such as with linear quadratic regulator (LQR) [6], Fuzzy-based PID [27], [28]. 
The PID was hard to be tuned to cope with the non-linear nature of the underwater environment, in which it 
typically produces high overshoot and high steady-state error. Due to the limitation of PID, artificial 
intelligent (AI) based controllers such as fuzzy logic controller (FLC) and artificial neural network (ANN) 
had been introduced to control ROV. 

ANN was used to predict the performance of the ROV depth system based on previous input and 
minimize the cost function [29]. Then, the best input is suggested. The ANN results were superior compared 
to other controllers that were experimented with. ANN was also used to tune the PID and adapt with the 
depth changing of ROV [30]. It was also used for the same purpose in [31] by implementing a radial basis 
function neural network (RBFNN) for trajectory tracking for the autonomous underwater vehicle (AUV). 
Both showed good results. However, the downside of ANN was the long computational time that may lead to 
a lagging problem. Another AlI-based controller for the ROV system is the FLC which was applied in ROV 
[21], [23], and AUV [32]. The FLC controller can cope with an unknown mathematical modeling system. 
Implementation of FLC eases the need for precise and complex hydrodynamic modeling of the vehicle. For 
example, FLC was successfully used to tune the PID controller for underwater vehicle [33]. 

However, even with the adaptability advantage, FLC has a challenging level of complexity. Hence, 
a simplified single input fuzzy logic controller (SIFLC) is proposed to control the depth of ROV. Previous 
works have shown that SIFLC has excellent performance, and it exactly resembled conventional FLC 
transient response [34], [35]. SIFLC reduces the input of conventional FLC into single input single output 
(SISO) system. Normally, a trial-and-error (heuristic) method was used to find the optimum parameter. 
Consequently, it takes more time execution to find the optimum parameters. 

This paper presents a new tuning approach of SIFLC with gradient descent algorithm (GDA) and 
particle swarm optimization (PSO) implementation for ROV depth control. The ROV was modeled using 
system identification to simulate the depth system. PID controller was applied to the model as a basic 
controller. SIFLC was then implemented in three different tuning approaches which were trial-and-error 
(heuristic), GDA, and PSO. The output transient was simulated using Matlab/Simulink and the percent 
overshoot (OS), rise time (Tr), and settling time (Ts) of the systems without and with controllers were 
compared and analyzed. In terms of depth control, the overshoot (%OS) is an important parameter to observe 
as a high value may damage the ROV or its investigation place [14], [18]-[20], [36]. The rise time (Tr) 
shows the time taken to get to the desired point while the settling time is the time ROV stabilizes at a steady 
state. 


2. SYSTEM MODELLING 

In this paper, the ROV was modeled using the system identification (SI) method in a Matlab 
computing environment. For system identification, the heave or vertical movement of ROV is being tested 
experimentally. Real-time input-output experimental data was gathered. 5 steps need to be considered in 
implementing system identification. Figure 1 shows the 5 steps for the SI approach. The steps are observation 
and data gathering, model structure selection, model estimation, model validation, and model application 
[37]. There are two sets of data needed for this experiment: training and validation data. The multi-sine signal 
was used to get the experimental data for training and validation. The input and output data were recorded, 
and the transfer function of the system was estimated using the Matlab toolbox. The input given to the ROV 
system can be a pulse, steps, random binary sequence (RBS), pseudo-random binary (PRBS), m-level 
pseudo-random (m-PRS), and multi-sine [37]. 

In the Matlab computing environment, multi-sine input was given to the system and the instrument 
variable (IV) approach was selected. Then, the selected model structure was implemented for model 
estimation and model validation to generate an ROV model. Finally, the model generated is used to design 
the ROV controller. The experiment was conducted in a controlled environment where the disturbances were 
not considered. In the instrument variable approach, the IV model was combined with 3 poles and 2 zeros for 
the transfer function. The best-fitting match was 96.43%. The transfer function generated is shown in (1). 
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The generated output transient response is shown in Figure 2. The output result had no overshoot, 
18.18 s of Tr, 33.21 s Ts, and 0.1947 of steady state error (SSE). Although the ROV did not overshoot, these 
results showed that the ROV model had a 19.47% target error and it took approximately half a minute to 
stabilize. This generated model was then simulated in Matlab/Simulink as a closed-loop system shown in the 
block diagram in Figure 3. 

This generated model was then simulated in Matlab/Simulink as a closed-loop system shown in the 
block diagram in Figure 3. From Figure 4, the closed-loop model had a faster Tr (9.07 s) and Ts (14.76 s) 
compared to the open-loop result but the steady-state error shoot up to 55.55% from the input given to the 
system. From the output result, the controller needs to be applied to get a better output response. 


Observation and data gathering 
Model structure selection 


Model estimation 


Model validation 


Model application 


Figure 1. System identification approach for modelling of ROV 


RiseTime: 18.1810 
SettlingTime: 33.2064 
SettlingMin: 0.7269 
SettlingMax: 0.8053 

Overshoot: 0 
Undershoot: 0 
Peak: 0.8053 
PeakTime: 71.6735 
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Figure 2. Transient response of the ROV model 
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Figure 3. Matlab/Simulink closed-loop block diagram 


Indonesian J Elec Eng & Comp Sci, Vol. 25, No. 3, March 2022: 1344-1355 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 1347 


oof + 7 7 + 7 + 7 T 


osp an 

0.7} fo 

f 

osp f 
© OS} Vá + ji t t 
2 dan Se SS SF 
= 04 f : t I i 
a 


OPEN LOOP SYSTEM 


| #-@@ CLOSE LOOP SYSTEM 


Time 


Figure 4. Open-loop and close loop transient output comparison 


3. CONTROLLER DESIGN 

In this paper, the PID controller was designed using auto-tuning provided by Matlab/Simulink. The 
SIFLC controller was designed and tuned using a heuristic method, GDA, and PSO. PID was used as a basic 
controller to be compared with the SIFLC controller designed. 


3.1. PID controller 

PID controller is the basic controller applied to the ROV system. The P, I, and D blocks were put in 
parallel in front of the plant to control the system. The P-block counters the direct error, the I-block rectifies 
the total errors in the system while the D-block minimizes the speed of the errors. The P-controller will make 
the response faster but intends to produce overshoot. The I-controller eliminates SSE while the D controller 
decrease overshoot. The PID controller block diagram is shown in Figure 5. The PID was tuned using 
automatic tuning in Matlab/Simulink [38]. 


a m (E) 


OPEN LOOP SYSTEM | 


PID Lae 
fo Lae 


Figure 5. PID controller block diagram 


3.2. FLC controller 

Fuzzy logic controller (FLC) is a human decision based controller introduced by Lotfi A. Zadeh. It 
was introduced in 1965. What do be done or reaction of the system is based on human perspective of the thig 
itself. In FLC algorithm, there are four (4) basic components; fuzzification, knowledge based (rules), 
inference engine and defuzzification. The fuzzification change the raw data into membership functioned, 
knowledge based create rules for decision making, inference engine which act as intelligent system and 
defuzzification change the fuzzy decision data into real output raw data. Figure 6 shows the basic 
configuration of FLC component in a block diagram [39]. 

The common 7 X 7 FLC table is shown in Table 1. The Z in the table stand for zero which indicate 
the center of the decision table. PL stand for positive large, PM stand for positive medium and PS stand for 
positive small. The NL, NM, and NS are the opposite of PL, PM and PS where they indicate the negative side 
of the table. The output result for FLC is selected based on the table. 


| 
based | 


Figure 6. Basic configuration of FLC component 
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Table 1. 7 X 7 FLC table 


Err vs du/dt or 1/s PL PM PS Z NS NM NL 


LNS 


3.3. SIFLC controller 

SIFLC controller is designed based on conventional FLC designed. The conventional FLC table, 
Table 1 is manipulated using signed distance method (SDM) which reduced the rules table to a one- 
dimensional array [40], [41]. From Table 1, it can be seen there is a consistent pattern in the decision-making 
of the FLC output. From Tables 1 and 2 diagonal lines were created which are named LZ and LNS. ‘d’ is the 
distance between LZ and LNS given by (2). The lambda (à) equation is shown in (3) and (4). 


_ wtZedA ow Zea 

~ i422 1422 t V1+22 (2) 
é+de=0 (3) 
> t (4) 


e 


Figure 7 shows the derivation of d, which is the distance between point, Q, and point, P. The 
conventional FLC table is now reduced to Table 2 where the diagonal line was represented by LNL, LNM, 
LNS, LZ, LPS, LPM, and LPL while NL, NM, NS, Z, PS, PM, and PL represent the output of corresponding 
diagonal lines. 


Main Diagonal line, L, 


INS 


é+Aae =O 


Figure 7. Derivation of d, the distance between point Q and P 


This input-output of SIFLC can be replaced by a lookup table. SIFLC was then tuned using the 
proposed lambda (4) tuning method. The value of (A) was heuristically changed to obtain the best output. The 
(A) was linked to the FLC by the input of the FLC. The range of error and integral error was plotted in a 
graph shown in Figure 8. 


Table 2. Reduced FLC table using SDM 
d LNL LNM LNS LZ LPS LPM LPL 
out NL NM NS Z PS PM PL 
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Figure 8. Plotted graph of input 2 versus input 1 FLC 


3.4. SIFLC heuristic tuning method 

The gradient of the line is lambda (A). The lambda (A) was varied up and down heuristically. The 
variation of (A)-tuned SIFLC result was then analyzed, and the best result was selected. Figure 9 shows the 
flow diagram of the heuristic tuning process. 


Increase and decrease the 
value of (A) linearly 


Record the (4) value 


Figure 9. Flow diagram for SIFLC heuristic tuning 


3.5. SIFLC GDA tuning method 

GDA is an algorithm that iteratively runs until it reaches the minimum value of a function. The 
GDA is used to replace the heuristic lambda (A) tuning for SIFLC. The objective function was obtained from 
the predicted output when compared to the input given. It is a simple mathematical method that is based on 
the differential equation where the initial point output was moved towards the targeted output by calculating 
the errors. Derivative of an objective function will determine the weight of the objective function for the next 
point. Two (2) important parameters are considered which are the direction of movement and the size of the 
step. The direction of movement is defined by the tangent of the initial point. The steepness of the tangent 
line also shows how near the point is to the minimum point and influences the choice of the learning rate. 
Figure 10 shows the flow diagram of the gradient descent algorithm [42]. The basic equation is shown in (5). 


P1 = P? — (Err * L) (5) 


Where, 

P? = current position 
Pt = next position 
Err = error (P? — dP?) 
L = Learning rate 
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Determine step size. 
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Figure 10. Flow diagram of gradient descent algorithm 


3.6. SIFLC PSO tuning method 

PSO was proposed by [43] in 1995. It is inspired by the behaviors of fishes schooling and birds 
flocking to search for foodstuff at a certain speed and position. The likeness is recognized between a particle 
and a swarm element [34], [44]. The particle movement is categorized by two factors: its current position x 
and velocity v, respectively. It has been useful effectively to a variety of optimization problems [45]-[47]. 
The particle swarm optimization algorithm is analyzed by using standard results from the dynamic theory 
[48]. The PSO algorithm begins by initializing the swarm randomly in the search space. Two consecutive 
iterations, t and t + 1 correspond to the position x of each particle that changes during the iterations by adding 
a new velocity v. The new velocity is estimated by summing an increment to the previous velocity value. The 
increment is a function of two components representing cognitive and social knowledge [49]. The cognitive 
knowledge of each particle is included by evaluating the difference between the current position x and its 
best position, pbest. The social knowledge of each particle is incorporated through the difference between its 
current position x and the best swarm global position achieved, gbest. He cognitive and social knowledge 
factors are multiplied by randomly uniform generated terms respectively [49]. Equation (6) shows the 
position vector while (7) shows the velocity vector. P in the equation is pbest while G is gbest. 


XPT = XP + 4 (6) 
Viti = wV? + arn (PÈ E xf) + C212 (G z xf) (7) 


Xt+1 = next position 
mae . 
Xt = current position 
a ; 

Vt+1 = velocity 


= 
wV =enertia (maintain current movement direction) 


w = weight (increase and decrease exploitation and exploration) 
Pt = personal best position 
Gt = group best position 


C,&Cz = the impact factors 
7, &rz = the random value (0 — 1) 
t = number of eteration 


4. RESULTS 

The output for all controllers designed was combined into one block diagram to compare the result. 
There are 6 signals analyzed which are step input, open-loop, closed-loop, PID, SIFLC heuristic, SIFLC 
GDA, and SIFLC PSO. Figure 11 shows the block diagram for the 6 signals investigated. From the block 
diagram, Scope 1 shows the 6 signals while scope 2 is used to compare between PSO result when the 
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equation was used (SIFLC PSO) and PSO result when the lookup table was used (SIFLC PSO1). This scope 
output shows an identical result (Figure 12) as expected. Figure 12 shows the output result for scope 1. In 
Figure 13, SIFLC GDA showed the most identical result to the step input given. It was then followed by the 
SIFLC heuristic. The SIFLC PSO showed improvement in the Tr but with a slight steady-state error. The PID 
showed some overshoot but no steady-state error. The output result is tabulated in Table 3. 
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CLOSE LOOP SYSTEM 
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Figure 11. Block diagram for the 6 signals investigated 
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Figure 12. Comparison result between PSO result using command windows (SIFLC PSO) and Simulink 
(SIFLC PSO1) 


New lambda tuning approach of single input fuzzy logic using gradient descent ... (Fauzal Naim Zohedi) 


1352. O ISSN: 2502-4752 


T 
To. 
E 
= 
> 
0. 
a Siep 
—@— OPEN LOOP SYSTEM 
02 CLOSE LOOP SYSTEM 


— — — PID CONTROLLER 
SIFLC GDA 


SIFLC Heuristic 
A— SIFLC PSO 


0 5 10 15 20 25 30 35 40 45 50 
Time (seconds) 


Figure 13. The output result of scope 1 


Table 3. Output result of the controllers' implementation to ROV system 
PID SIFLC Heuristic _SIFLC GDA SIFLC PSO 


Tr(s) 7.0665 7.2529 0.7992 2.3686 

Ts (s) 24.6687 10.9736 0.9790 12.2348 
%OS 7.3613 0.7988 0.1021 16.2368 
SSE 0 0 0 0.1 


From the bar chart in Figure 14, it is obvious that the SIFLC GDA control method showed the best 
result as it managed to get the lowest values for all the performance parameters. In terms of the rise time, Tr, 
the SIFLC GDA recorded 0.7992 s while the SIFLC heuristic was the worst at 7.2592 s, nearly ten times the 
SIFLC GDA record. The heuristic approach and PID approach had almost similar values at 7 s. These 
records showed that SIFLC GDA had the fastest response time to the change of input level. For settling time, 
Ts, the SIFLC GDA method recorded 0.9790 s and PID recorded the worst at 24.6687 s. Other than SIFLC 
GDA, other methods to stabilize the ROV took ten to twenty times longer. For the percentage of overshoot, 
the SIFLC GDA showed the best result which is at 0.1021%. It is then followed by the SIFLC heuristic 
(0.7988%), PID (7.3613%), and SIFLC PSO (16.2368%). These results showed that when the ROV was 
controlled using the SIFLC GDA method, the error to attain the target was very small. Figure 15 shows the 
comparison of SIFLC GDA with step signal. The SIFLC GDA looks nearly identical to the given step input. 
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Figure 14. Bar chart of the performance parameters 
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Figure 15. Comparison of SIFLC GDA with step signal 


5. CONCLUSION 

Two controllers with different tuning had been applied to the ROV depth system. A new tuning 
approach of SIFLC controller based on lambda (A) is proposed and was compared with the basic PID 
controller. The SIFLC GDA showed the best result as it had the lowest values in all performance parameters 
investigated which were 0.1021% (OS), 0.7992 s (Tr) and 0.9790 s (Ts). The SIFLC PSO suffered from high 
overshoot and some steady-state error. The basic PID controller also suffered from some overshoot and had a 
long settling time. The SIFLC Heuristic had a better result compared to the PID controller in Ts and %OS 
however it was difficult to tune and required experience and time. The SIFLC GDA was able to obtain 
plausible results because it uses specific tuning of the objective function which is based on all parameters; 
%OS, Tr and Ts. On the other hand, the SIFLC PSO had a larger error compared to SIFLC GDA because the 
objective function used was the absolute mean error value. It doesn’t specifically tune based on each 
parameter. From all results, it is proven that SIFLC lambda (A) tuning approach successfully produced good 
output result. With the implementation of an optimization approach such as GDA and PSO, a better output 
result can be obtained. This happen because it optimizes the value of lambda (A) to suit the input output 
result. The objective function selected in running the optimization approach also plays an important role in 
getting a good result. For future implementation, varieties of objective functions implementation can be 
studied and proposed to the system. 
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